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We consider a class of late-decaying dark-matter models, in which a dark matter particle decays 
to a heavy stable daughter of approximately the same mass, together with one or more relativistic 
particles which carry away only a small fraction of the parent rest mass. Such decays can affect 
galactic halo structure and evolution, and have been invoked as a remedy to some of the small-scale 
structure-formation problems of cold dark matter. There are existing stringent limits on the dark 
matter lifetime if the decays produce photons. By considering examples in which the relativistic 
decay products instead consist of neutrinos or electron-position pairs, we derive stringent limits on 
these scenarios for a wide range of dark matter masses. We thus eliminate a sizable portion of the 
parameter space for these late-decay models if the dominant decay channel involves Standard Model 
final states. 

PACS numbers: 95.35.+d, 95.85.Pw, 98.62.Gq, 98.70.Vc 



I. INTRODUCTION 

Despite the ample gravitational evidence for the exis- 
tence of dark matter (DM), its nature remains unknown. 
A promising strategy for identifying the dark-matter par- 
ticle is to search for a signature produced by DM decay 
or annihilation, inside the Galactic halo or at cosmolog- 
ical distances. Although such a signature has not yet 
been identified, the observed radiation backgrounds have 
served to constrain various decay and annihilation chan- 
nels [j|4la|. In many cases, this has resulted in stringent 
limits on specific DM candidates. 

In this paper we derive constraints for a class of models 
in which DM decays dominantly via the process 



X 



x' + i, 



(1) 



where \' is a massive stable particle, and I denotes one 
or more light (or possibly massless) particles. Such decay 
modes have been discussed in Refs. |14l - ll9l |. We focus on 
the case where the mass splitting between the unstable 
parent \ and the daughter particle \' is small, 



Am 



£TOv 



(2) 



where e -C 1. Since the decay replaces an x DM particle 
with an \' DM particle of approximately the same mass, 
it does not change the DM energy density significantly 
and it is thus not constrained to occur at time scales much 
larger than the age of the Universe. This is quite different 
from decaying DM scenarios in which the entire energy 
of the DM is converted into relativistic species during 
the decay, such as the model of Ref. [20[ . The latter sce- 
nario is stringently constrained by the cosmic microwave 
background anisotropics, which set a lower limit on DM 
lifetime of 123 Gyr at la confidence level HJ. The DM 
decay models we consider in this paper, in which most 
of the energy is retained in the form of nonrelativistic 
matter, may be most effectively constrained by compar- 
ing the flux of the light particle(s) I emitted in the decay 
against the observed radiation backgrounds. 
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FIG. 1: Constraints on the dark- matter decay channels: (i) 
X — > X + 7) using the 7-ray line emission limits from the 
Galactic center region and the isotropic diffuse photon back- 
ground (irregular and smooth solid lines, respectively) (ii) 
X — > x' + e ~ + e+ (dashed line); (iii) x x' + v + v (dotted 
line). The regions below the lines are excluded. 



An interesting feature of this class of models is that 
they offer a possible solution to some of the discrepancies 
between galactic observations and small-scale structure 
predictions of cold dark matter (CDM) simulations. In 
CDM simulations, clustering proceeds hierarchically, in a 
"bottom-up" fashion, resulting in rich structure at small 
scales. Several disparities with galactic observations have 
been identified |22h35| . For the disagreement to be re- 
solved within the CDM paradigm, some mechanism that 
modifies the standard CDM structure-formation picture 
is necessary. If dark matter decays according to Eq. ([T]) 
at time scales comparable to the age of the Universe, the 
kinetic energy acquired by the heavy daughter, will ef- 
fectively heat up the dark-matter haloes, and cause them 
to expand. This can alleviate two of the most glaring 
CDM problems, the cuspy density profiles of dark-matter 
haloes and the over-prediction of satellite galaxies [T3 - 
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m. 

Yukscl and Kistler [1] obtained bounds for the decay 
channel x ~> x' + 7; using observations by SPI, COMP- 
TEL and EGRET. Since photons are the most easily 
detectable particles of the Standard Model (SM), these 
bounds constrain late-decaying DM scenarios in the most 
stringent fashion. However, the branching ratio for direct 
photon production may be tiny in many models. In this 
paper we investigate DM decay modes to other possible 
SM particles. We focus on two cases: when the light 
daughter particles I produced in the decay are a pair of 
the lightest charged leptons of the SM 

X^X' + e~ +e+ (3) 

or a pair of the least detectable stable SM particles 

X -> x' + v + v. (4) 

Given that photons provide the most stringent bounds, 
and neutrinos the least stringent, our results, together 
with those of span the full range of constraints for 
all SM final states. Moreover, constraining the decay 
into neutrinos places a general lower bound on the DM 
lifetime for decay to any SM particle. 

The production of electrons and positrons by dark- 
matter decay can be constrained both by the observed 
positron flux and by the observed photon background. 
Electrons and positrons propagate in the Galaxy and 
produce photons in a variety of ways. We consider a 
comprehensive list of effects which result in photon emis- 
sion associated with the dark-matter decay into pairs, 
including inverse Compton scattering, bremsstrahlung, 
synchrotron radiation, in-flight and at-rest annihilations 
and internal bremsstrahlung. We compare the expected 
signals with data from Fermi LAT, COMPTEL, EGRET 
and INTEGRAL. The decay into vv can be constrained 
by the atmospheric-neutrino background, and by data 
from the Super-Kamiokande search for the diffuse super- 
nova neutrino background. 

In Sec. [II] we derive the constraints on the decay chan- 
nels ([3]) and (UJ). We present our results in Fig. [1] to- 
gether with the constraints derived in Ref. [l| for decay 
to a photon. In Sec. IIII1 we investigate whether these 
limits permit DM decay to have a significant effect on 
the formation of galactic structure. We focus on late- 
decay scenarios, at times > 0.1 Gyr, in accordance with 
the work of Refs. [lB - [l9l ]. For the decay modes described 
in Eq. ([3]) , and Eq. ((4]) , our constraints eliminate a sizable 
band of the r vs m x parameter space. However, signifi- 
cant parameter space for which DM decay may influence 
halo structure remains available. In Sec. IIVI we turn our 
attention to models which are not necessarily motivated 
by the structure formation debate, nevertheless they in- 
volve decay of a relic population of particles into states 
of similar mass. In such models, the decaying particles 
may amount only to a small fraction of the total DM den- 
sity of the Universe. We explore the applicability of our 
constraints on these scenarios and comment on specific 
particle-physics models. 



II. CONSTRAINTS ON DARK-MATTER 
DECAY 

The energy of the light particles produced in the decays 
depends on the mass splitting between the parent and 
the heavy daughter particle, Am. Their flux is inversely 
proportional to the lifetime, r, and the mass, m x , of 
the parent particle. The latter determines the number 
density of the DM particles. Observations thus set a 
lower limit on the product m x T as a function of Am. 

However, if r is comparable to the age of the Uni- 
verse To ~ 13 Gyr, an essential element in the pro- 
posed late-decaying DM resolutions of the small-scale 
structure problems of CDM [TfHH, 

the abundance of 

X particles may have changed significantly due to decay. 
Observational bounds should then be cast in the form 
m x T e TQ / T ^ J(Am), reflecting the exponential suppres- 
sion of the parent-particle number density today. (The 
observational bound, J"(Am), is presented in Fig. Q] and 
will be discussed further below.) Note that for sufficiently 
short lifetime, t « To, all the x particles decay to the 
stable x' state very early on, thus leaving no observa- 
tional signal in the Universe today, while for sufficiently 
long lifetimes, r > ro, the current decay probability will 
be very small. The lifetimes of interest to us, i.e., those 
that we can probe via indirect detection of the DM decay 
products, lie between these two extremes. 

In our analysis below, we will assume that the x P ar ~ 
tides originally constitute all of the DM, with no signif- 
icant admixture of the similar-mass x' particles. Com- 
parable initial abundances of x an d x' might be antic- 
ipated in canonical thermal-relic realizations of such a 
DM model. However, our limits can easily be rescaled to 
account for multicomponent DM, and a down-scaling of 
our constraints by a factor of ~ 2 would not qualitatively 
affect our conclusions. 

The massive daughter particle x' wm be produced non- 
relativistically, provided that Am -C m x . Then, to a 
good approximation, the light particles share the energy 
Am released in the decay, such that E\ + E% ~ Am, 
where E\,E<z are the energies of the two light leptons. 
For the 3-body decay channels described in Eq.® and 
Eq.((4]), the energy spectra of the daughter particles is in 
general model-dependent. However, we find that the con- 
straints on the dark-matter decay are rather insensitive 
to the actual spectrum. We illustrate this by present- 
ing results for two limiting cases of the daughter energy 
spectra: 

i. Monoenergetic, where the energy spectrum per de- 
cay for each light daughter species is 

dN/dE = S(E - Am/2). (5) 

ii. Flat, where is spectrum is uniform over the allowed 
energy range such that 

dN/dE = l/(Am - 2m/), (6) 
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where mi = m e for the decay into e and mi — 
for the decay into neutrinos. 

For each of the decay modes analyzed, we pick the most 
conservative limits to further continue our discussion 1 . 

The constraints on DM decay depend on the DM halo 
density profile. Although there are uncertainties in the 
profile, particularly near the Galactic center (GC), we 
shall see below that these uncertainties are subdominant. 
For consistency with the results of other analyses which 
we shall employ in our calculations, we use the Navarro- 
Frenk- White (NFW) [H[ profile throughout. The stan- 
dard parametrization is 



Iy > (r/r.)T[l + (r/r,)°'](/'-7)/«' w 

where (a, 7) = (1,3,1). For the Milky Way r s — 
20 kpc, and the normalization po is fixed such that at 
the solar-circle distance J? sc = 8.5 kpc, the density is 
p(R BC ) = 0.3GeVcm- 3 . 



A. Decay into electrons and positrons. 




Am [GeV] 

FIG. 2: Limits on dark-matter decay process \ X + 
e~ + e + , obtained from internal bremsstrahlung (IB, cyan), 
in-flight annihilation (IA, blue), annihilation at rest (511 
keV, purple), inverse Compton scattering and bremsstrahlung 
emission (IC, green), synchrotron radiation (Synch, orange), 
and positron flux (e + flux, red). Solid lines correspond to 
a monoenergetic spectrum of electrons and positrons, while 
dashed lines assume a flat injection spectrum over the allowed 
energy range. The regions below the lines are excluded. 



Electrons and positrons injected in the Galaxy will pro- 
duce photon fluxes spanning a wide range of energies, by 
a variety of different mechanisms. The injection rate can 
thus be constrained by the resulting contribution to the 
Galactic photon backgrounds. We will also derive con- 
straints based upon direct detection of positrons. Below, 
we describe the methods used to estimate the signal pro- 
duced by each effect. We present the resulting constraints 
in Fig. dJ 



uncertainties associated with the astrophysical environ- 
ment. IB is thus a very clean and reliable technique for 
obtaining limits on DM decay to charged particles, with 
the only source of uncertainty residing in the DAI profile 
adopted. 

The IB differential decay width is [37} 



dT 



IB 



dE 



= T x h(E), 



(8) 



1. Internal bremsstrahlung. 

Electromagnetic radiative corrections to the DM de- 
cay channel to charged particles inevitably produce real 
7 rays. If the process x ~ ► x' + e± takes place, then 
X — > x' + e± + 7 must also take place since either one 
of the final state charged particles can radiate a photon. 
This process is called internal bremsstrahlung (IB). Note 
that the photon in the IB process arises directly in the 
Feynman diagram for the decay, and is unrelated to reg- 
ular bremsstrahlung which results from interaction in a 
medium. Importantly, the IB spectrum and normaliza- 
tion is approximately model-independent. In addition, 
as the IB photon flux does not depend on the propaga- 
tion of charged decay products in the Galaxy, it is free of 



where E is the energy of the photon emitted, T = t _1 
is the rate of lowest order decay process x —> x' + e± i 
and h(E) is the photon spectrum per x ~> x' + e± decay. 
This spectrum is independent of the new physics which 
mediates the decay, and is given by 



^)^Iln(4 
7r E \mi 



1 + 1 7 



(9) 



where s' = s(l — E/E max ), with s being the energy 
released in the decay, and E max the maximum energy 
that can be imparted to the radiated photon 2 . In the 
3-body decay we consider here, s = (Am — 2m e ) 2 , and 
£max = - e/2) ~ ^ 3 



As a consistency check, we verify that for a specific Am, the limit 
obtained on m x Te T o/' T when an extended injection distribution is 
used, does not exceed the maximum value of this limit obtained 
for monoenergetic injection for mass splittings up to 2Am, i.e. 
[m x reWT( Am )] ^ ma Xfl<2Am [m x r e Ta ' T {n)] 



2 These general definitions of s, s' reproduce the correct expres- 
sions for the case of annihilation into two charged particles [37^1 
and for muon decay [38ll . 

3 Note that the phase space available to the photon when one of the 
decay products is emitted nonrelativistically (i.e. when e < 1) 
is larger than when all decay products are relativistic (e.g. muon 
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The shape of the IB photon spectrum is determined 
by the 1/E factor in Eq. ©, while the logarithmic factor 
provides a sharp cutoff at E ~ E max . Since the galac- 
tic photon background varies approximately as E~ 2 , the 
high-energy region of the IB photon spectrum will pro- 
vide the most stringent constraints. 

The flux of photons due to IB will be proportional to 
the dark matter density integrated over the line of sight. 
At an angle ip from the Galactic center, the line-of-sight 
integral is 



dl p x - 21R SC cos?/>-M 2 ) , (10) 



where / max = y-R^w _ ^sc sin i 3 + -^sc sin?/'. The aver- 
age of J over a cone of half-angle tp centered on the GC 
is 

J ^=MlJ Q WsinW, (11) 

where Af2 = 2n(l — costp). 

The expected photon flux, per unit solid angle, due to 
IB is then 



1 



dE Air m x Te T °/ T 



J±n HE). 



(12) 



COMPTEL observations show that the 7-ray back- 
ground remains approximately constant within |^| < 30°. 
Thus, along the Galactic plane, the signal-to-background 
ratio for IB is maximized at the GC. Here, we will average 
over the region tp < 5°, which exceeds the angular res- 
olution with which both COMPTEL and EGRET data 
have been reported [H, H(| ■ The IB limits are shown in 
Fig. [2] (cyan line). 



2. In-flight and at-rest annihilations. 

Positrons produced by DM decay will propagate in the 
interstellar medium, lose energy, and annihilate with elec- 
trons. Most positrons will survive until they become non- 
relativistic, and annihilate at rest to produce a 511 keV 
annihilation line (together with a three-photon contin- 
uum). However, a portion of the positrons annihilate in 
flight while still relativistic. These in-flight annihilations 
(IA) can produce a significant flux of 7 rays [4l[ . 

Reference [4l[ considered the injection of positrons at 
the Galactic center, and related the flux of IA 7-rays to 
the flux of 511 keV photons. This was used to derive a 



decay, where e ~ 1). This is because a heavy daughter particle 
can absorb recoil momentum, without carrying away a significant 
amount of energy. 



constraint on the injection energy of positrons that con- 
tribute to the observed Galactic 511 keV line. 

We will follow a similar approach to obtain limits on 
DM decay into e . In our case, however, the overall nor- 
malization of the photon flux is not fixed by the observed 
511 keV line intensity, but is instead determined by the 
DM decay rate. The continuum 7-ray flux due to IA, and 
the 511 keV photons from annihilation at rest, thus yield 
thus two separate sets of constraints. The positron in- 
jection, deceleration and annihilation budget is described 
by the equation [42[ 



dn d ( d"f 



dt <9 7 V dt 



-n H cr(7)/3n(7,t) 



dN , s 

r inj -, (13) 



where n(^,t)d^ is the density of positrons with energy 
7 m e , at time t, T^j = p x e~ T °/ T /m x T is the positron- 
density injection rate due to DM decay, and dN/dj is the 
positron spectrum per decay. The positron annihilation 
cross section 17(7) on electrons at rest is 43] 



7+1 



7 2 + 4 7 + 1 



In 



1 



,--(7 + V^l)-#= 

7-1 v > x/l 2 - 1 

(14) 

where r e is the classical electron radius. 

At energies lower than ~ 1 GeV, ionization and 
Coulomb losses (in a neutral or an ionized medium, re- 
spectively) dominate over energy loss due to synchrotron, 
bremsstrahlung and inverse Compton effects, while at 
higher energies synchrotron losses are more important. 
The energy loss rates due to ionization, Coulomb scat- 
tering and synchrotron emission are [42| 



d" 

G?7 

~dl 

d~i 
~dt 



4.4- 10- 15 s- 1 f—!^- 

V0.1 cirr 

1.5- lO-^s- 1 (—^ 

\0.1 cn 



cm 
U B 



ln7 + 6.8 
ln 7 -r-75.9 



(15) 



(16) 



0.2 



eV 



(7 2 "1), (17) 



where n a is the hydrogen number density of the medium, 
and Ub = B 2 /8w is the energy density of the magnetic 
field. Equations (ITS"]) and Eq. (fT6)) apply to a fully neutral 
and or fully ionized medium, respectively. As a conser- 
vative choice, we set the ionization fraction of the inter- 
stellar gas to be Xi 0.51. Compared with a completely 
neutral medium, energy losses in an ionized medium are 
larger and thus IA constraints weaker, by a factor of a 
few [13 

As evident from Eqs. (|15l) and it takes approxi- 

mately 10 7 yr for relativistic positrons to slow down and 
annihilate. If positron injection (DM decay) has been 
taking place over a much larger time period, t > 0.1 Gyr, 
the positrons reach a steady-state distribution. The time- 
independent solution of Eq. (ri"3]) is 



71(7) 



Px 



1 



m Y T e 



t /t 



1 , , dN „ 



(18) 
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where 



7 >7 



exp 





7")/3" 




d-y" 






dt 





dj" 



(19) 



is the survival probability, introduced in |4l| . of positrons 
injected at energy 7'm e to have not annihilated by the 
time they reach energy -ym e . 

Here, we have ignored the effects of spatial diffusion of 
positrons, which would require detailed modeling. This 
approximation is sufficient provided that the positrons 
remain trapped within the region we consider 4 . In what 
follows, we will calculate the expected photon flux due to 
annihilation in a cylindrical region of radius r max = 3kpc 
and half-height z max = 0.5 kpc, centered at the GC. This 
is well inside what is taken to be the diffusion area in most 
models simulating cosmic-ray propagation in the Galaxy. 
Because of the Galactic magnetic field, charged particles 
are largely expected to follow the magnetic lines and con- 
verge toward the GC, as long as their Larmor radius, 
rL = -ym/qB, is much smaller than the scale of variation 
of the magnetic field. In the Milky Way, B > few /iG, 
and (din B j dr)^ 1 ~ few kpc [5j. Positrons are thus ex- 
pected to follow the magnetic lines, even for energies 
much higher than the ones considered in this analysis. 
Any influx of positrons in the region under consideration 
would only lead to more stringent constraints than the 
ones derived here. (The possibility of transport of the 
positrons from the disk to the GC has been discussed in 
detail in 0].) 

The diffuse 7-ray background of the Galactic plane 
has been measured by COMPTEL in the energy range 
1-30 MeV for the region \l\ < 30°, |6| < 5° (Galactic 
longitude and latitude, respectively) and by EGRET, 
in the e nerg y range 30 MeV-30 GeV and the region 
|b| < 10° [3!J,|4(J. The cylindrical volume defined above is 
covered by these observations, and encompassed within 
a conical patch of solid angle Aft = 4i max sin b max ~ 
0.13 sr, with Z max = arcsin(r max /i? sc ) ~ 20° and 6 max = 
arctan[z max /(i? sc — r max )] ~ 5°. The comparison of the 
signal from the cylindrical region with the background ra- 
diation emanating from the larger conical region, keeps 
our constraints conservative. 

The flux of photons (per unit solid angle) due to IA is 
d$iA 1 1 f „ 7 f , / \ da{k,i) 



(Ik 



I I f f < 

M}^Rj c J dV J^ n ^<- 



dk 



4 Along with spatial diffusion, we also ignore the related effect 
of reacceleration. This becomes more important with increas- 
ing energy, since the diffusion coefficient in momentum space is 
D pp oc p 2 ~ s with S < 1 @, IEE3l- However, the probability 
to annihilate drops sharply at high energies, and thus pumping 
more energy in already energetic positrons does not affect the IA 
constraints significantly. 



where the spatial integration is over the cylindrical region 
considered. The differential cross section dCT ^,' 7 - ) for a 
positron of energy -ym e to produce a photon of energy 
E = k 



m„ is 



HE 



da(k 1 j) irrl 



dk 



7 2 - 1 



3+7 , 3+7 _ 1 

' 1+7 ' k W 



1 



k 
1+7 



(20) 



while the minimum positron energy 7 m i n m e required to 
produce a photon of energy km e is given by 



1 



7min(fc) — k — 



1 



2 2(2fc-l) 



(21) 



For the steady-state positron distribution 71(7), given in 
Eq. (fTB)) . the expected 7-ray signal from IA is 



IA 



dk 
where 

Q(fc,7) = n H 



1 fdV Px 
m x r e r °/ T 47ri?2 c Af2 



C 00 r} AT 

/ d7^Q(^,7), 

(22) 



7 v^i^y)^ (23) 



d-y 

7m ln (fc) 



7 



dk 



or 



The IA flux of Eq. ([22]) is compared to the 7-ray back- 
ground reported by COMPTEL and EGRET [Hi Ejj, 
and compiled in Ref. [l2[ • The resulting lower limit on 
m x T e r °/ T is presented in Fig. [2] (blue line) and increases 
with increasing Am, since positrons injected at higher en- 
ergy spend more time as relativistic, and are more likely 
to annihilate in flight. However, the annihilation prob- 
ability diminishes at high energies, which is manifest by 
the plateau in the lifetime bound for energies > 3 GeV. 

The majority of positrons survive until they become 
nonrelativistic. Then, they either annihilate directly with 
electrons to produce 511 keV photons, or they form a 
positronium bound state. Positronium subsequently an- 
nihilates, with probability 25% into two 511 keV pho- 
tons, or with probability 75% into a three photon contin- 
uum. The positronium fraction in the Galaxy is fixed 
by the relative intensities of the 7-ray continuum be- 
low 511 keV, and the 511 keV line. In our Galaxy 
/ = 0.967 ±0.022 [13. 

If DM decays into pairs, the flux of 511 keV pho- 
tons produced is determined, in the steady-state regime, 
by the rate at which positrons arrive at rest. For each 
nonrelativistic positron, there will be 2(1 — /) + 2//4 = 
2(1 — 3//4) photons contributing to the line. The ex- 
pected flux of 511 keV photons is 



$511 = 



2(l-3//4) JdV Px 



<re ro/r Aft J 1 



' dN „ 



(24) 



The 7-ray observations by INTEGRAL have re- 
vealed a 511 keV line emission of intensity 0.94 x 
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10~ 3 phcm -2 s _1 within an integration region \l\ < 20°, 
|b| < 5° [11]. The average flux per steradian is then 
0.0077ph cm -2 s _1 sr -1 . The limits on DM decay arising 
from annihilation at rest are presented in Fig. [5] (purple 
line). They exhibit a very slight negative slope towards 
increasing Am, which accounts for the (small) portion 
of positrons annihilating in-flight when injected at high 
energies. 

The IA constraints are stronger than those for anni- 
hilation at rest for average positron injection energies 
(E) = Am/2 > 80 MeV. This is considerably higher 
than the corresponding value of ~ 3 MeV found by an 
appropriate analysis in Ref . [41| . The latter limit was de- 
termined by requiring that the IA 7-rays do not exceed 
30% of the observed background, instead of the more 
conservative 100% adopted here. However, the more im- 
portant factor in reconciling these two analyses is the dif- 
ferent signal regions considered. The analysis in Ref. [4l[ 
considered the conical region within a ~ 0.37 kpc radius 
of the GC, which contains the peak of the 511 keV an- 
nihilation signal. This is to be compared to our larger 
cylindrical volume (r max = 3 kpc, z max = 0.5 kpc) out- 
lined above. The effect of choosing the larger volume 
is to make the annihilation-at-rest constraints stronger, 
with respect to those from IA, than for the choice of a 
smaller region at the Galactic center. The Galactic cen- 
ter is not the optimal observation region to use in setting 
an annihilation-at-rest limit, because that is where the 
observed 511 keV background is highest. In addition to 
resulting in more sensitive annihilation-at-rest limits, the 
larger volume also renders us insensitive to possible ef- 
fects of diffusion. 

Notice also that the photon fluxes, Eqs. ([22]) and (|24|). 
are independent of the medium number density , as 
long as ionization and Coulomb losses dominate, for 
Am < few GeV. Since the probability to annihilate at 
even higher energies is insignificant, the results are rather 
insensitive to n H , as well as Ub, throughout the energy 
spectrum. 



3. Inverse Compton scattering, bremsstrahlung emission, 
synchrotron radiation. 

At high energies IA become rather rare, and electrons 
and positrons propagating in the Galaxy produce radia- 
tion more efficiently via other mechanisms. Gamma rays 
are produced by bremsstrahlung and inverse Compton 
(IC) scattering of low-energy photons, with the latter 
effect yielding the dominant contribution. Synchrotron 
emission in the Galactic magnetic field gives rise to a 
radio wave signal. 

Zhang et al. H, @ calculated the photon spectrum 
expected if electrons and positrons are injected in the 
Galaxy by DM decay. They assumed a monoenergetic 
injection spectrum, and modeled the propagation of e 
in the interstellar medium in detail, including the effects 
of spatial diffusion, convection, energy loss and reaccel- 



eration. They compared the resulting spectrum at Earth 
with observed backgrounds, and encoded their calcula- 
tions in response functions which can be utilized to ob- 
tain constraints when convoluted with the DM decay 
spectrum of a particular model. Here, we use the re- 
sponse functions derived in Ref. Q for synchrotron radi- 
ation, and in Ref. @ for IC and bremsstrahlung emission, 
to obtain constraints for the x ~ * x' + e± decay process. 

The main source of uncertainties in deriving the re- 
sponse functions is poor knowledge of the various as- 
trophysical parameters which determine the propagation 
of e in the Galaxy. The most significant contribution 
comes from the height of the diffusion zone. The Galactic 
magnetic field, although also quite poorly known, con- 
tributes subdominantly to the total uncertainty. This is 
true even for the synchrotron response functions, since al- 
though the magnetic field may be rather uncertain in the 
Galactic bulge, the directions which optimize the signal- 
to-background ratio point away from the GC. Different 
choices for the DM halo profile have a small influence 
on the response functions, since the flux is only propor- 
tional to the density (as opposed to the p 2 dependence 
applicable for DM annihilation) [|| Q . 

The IC and bremsstrahlung emission is constrained by 
using the Fermi LAT 7— ray maps derived in [4^ |. for 
the energy range 0.5 GeV-300 GeV. Because of reaccel- 
eration, this permits limits for e injection energies in a 
much wider range, 0.01 GeV — 10 4 GeV. The strength of 
the limits, of course, sharply diminishes in the low-energy 
part of this range. At energies above 50 GeV, the Fermi 
data may suffer from significant background contamina- 
tion, while point sources have not been subtracted. These 
effects lead to more conservative constraints. There are 
also various astrophysical contributions to the 7-ray flux, 
such as nucleus-nucleus photoproduction via 7r° decay, 
and IC and bremsstrahlung radiation from cosmic-ray 
electrons and positrons. In the spirit of setting conserva- 
tive limits on DM decay, we will use response functions 
where these foregrounds have not been subtracted. The 
patch of sky chosen to compare the expected signal to 
observations is \l\ < 20° and -18° < b < -10° @. 

The synchrotron response functions derived in Ref. [j| 
correspond to 408 MHz, 1.42 GHz and 23 GHz sky maps, 
and yield constraints for e ± injection energies in the 
range, 0.1 GeV — 10 4 GeV. The directions that optimize 
the signal-to-background ratio are different for each of 
the sky maps used, and are located close to, but not at, 
the GC. (Similar sensitivity may be achieved with the 
isotropic diffuse flux measured by Fermi; see Ref.[5(| for 
relevant limits on the process x — >■ e^.) 

The response functions reported all assume a NFW 
halo profile. The expected signal from particles pro- 
duced by DM decay is compared to observed photon flux 
including a 2a error. The dependence on the propaga- 
tion model, for both the synchrotron and the IC plus 
bremsstrahlung emission, becomes more significant at e ± 
injection energies lower than about 10 GeV. This is be- 
cause the limits in that range rely solely on the effect 
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of reacceleration, which depends strongly on the model 
adopted, while no actual low-energy background data 
are used. For the synchrotron constraints, we adopt the 
"DR" model @, which yields moderately conservative 
constraints along the entire energy range. For the IC plus 
bremsstrahlung constraints we adopt the "LI" model, 
whose response functions are most comprehensively re- 
ported in [6J. 

The limits on DM decay from IC plus bremsstrahlung 
radiation and synchrotron emission are shown in Fig. [5] 
(green and orange lines, respectively). 



4- Positron flux. 

The DM decay into pairs can be constrained by 
direct observations of the positron flux on Earth. Uti- 
lizing the same propagation models as the ones used for 
the synchrotron emission, Zhang et al. [j| constructed re- 
sponse functions, c omp aring the expected positron flux 
to PAMELA data [51(. The positron fraction reported 
by PAMELA was converted into positron flux using the 
e ± total flux observed by the Fermi telescope [52j . 

The positron-flux response functions are based on ob- 
servations of positrons in seven energy bands in the range 
10.17 GeV — 82.55 GeV. Depending on the model used, 
they can provide meaningful constraints for positron in- 
jection energies within the range 0.1 GeV— lOTeV. These 
constraints appear to prevail over the IB, IA, IC and 
synchrotron constraints at the high-energy part of the 
spectrum. In the same spirit as before, we adopt the 
DR model, and obtain fairly conservative constraints for 
the entire energy spectrum. The results are presented in 
Fig. |3 (red line). 

Much like the photon-emission response functions, 
the positron-flux response functions are sensitive to the 
adopted propagation model, which determines the diffu- 
sion and reacceleration of the injected positrons. How- 
ever, they are insensitive to the DM halo profile, since 
most of the positrons come from the local region, ~ 1 kpc 
from the Sun [I] 5 . 



B. Decay into neutrinos. 

Neutrinos are the least detectable stable SM parti- 
cles, hence constraints on DM decay (or annihilation) 
to neutrinos can be used to set conservative but ro- 
bust lower limits on the DM lifetime (or upper limit on 
the annihilation cross section) to any SM final state Q- 
Palomarez-Ruiz 0, Q obtained limits on the decay chan- 
nel x ~~ * v+v. Here we will adapt the analysis of Refs. 0- 
4], to constrain the decay channel x ~^ x' + v + This 




10 1 10 2 
Am [GeV] 



FIG. 3: Limits on dark-matter decay \ — > \ + v + v. As in 
Fig. [2] the solid (dashed) lines correspond to monoenergetic 
(flat) injection distribution of v, v. The regions below the 
lines are excluded. 



will set the most conservative limit on DM decay modes 
of the form x x' + h where / is any SM final state. 

The neutrino flux on Earth from DM decay in the 
Galactic halo is 



1 



1 



= 7 2dV 

dE u 4tt m y re T °/ T 3 dE v 



(25) 



5 This is consistent with our assumption in Sec. Ill A 21 that the 
do not escape from the inner part of the galaxy. 



The multiplicity of 2 accounts for the sum of v and v 
produced in the decays. As in Refs. 0, 0|, we assume 
equal decay width to all three neutrino flavors, which 
accounts for the factor of 1/3 in Eq. (|25|) . This is a 
reasonable approximation, since neutrino flavor oscilla- 
tions between the production and detection points will 
considerably weaken any preference toward a particular 
flavor. For simplicity, and in order to minimize the un- 
certainty arising from a particular choice of halo profile, 
we average the expected neutrino signal over the whole 
sky, ip = 180°. Directional information, whenever avail- 
able, is in general expected to lead to more stringent 
limits @,0,|1]. 

For energies E v > 50 MeV, the neutrino flux on 
Earth is dominated by atmospheric neutrinos. The at- 
mospheric neutrino flux has been well measured by a 
number of experiments [53l - [58j and is in good agree- 
ment with theoretical predictions. We utilize the results 
of the FLUKA [H^, [60] calculation of the atmospheric 
+ background, over the energy range 50 MeV to 
10 TeV. We determine a DM decay limit by requiring 
the expected neutrino signal from DM decay not to ex- 
ceed the atmospheric neutrino flux, integrated over en- 
ergy bins of width A log 10 E v ~ 0.3. This choice of bin 
size is in accordance with that adopted in Refs. [2|, |3j, 
and encompasses the experimental resolution of the neu- 
trino detectors. (We note that future measurement by 
IceCube+DeepCore will lead to sensitive limits at high 
energy [6l|, |62[.) 

At low energies E v < 100 MeV, the relevant data 
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come from the diffuse supernova neutrino background 
(DSNB) search performed by the SK experiment |63j . 
SK searched for positrons produced by incident v e on 
free protons inside the detector, via the inverse /3-decay 
reaction v e +p — > e + +n. Incoming v e and v e interact also 
with bound nucleons, producing electrons and positrons. 
The main sources of background for these observations 
are atmospheric v e and v e , and the Michel electrons and 
positrons from decays of subthreshold muons. A DSNB 
signal was not detected. In Refs. a x 2 analysis of 
these data was performed, and used to place limits, at 
90% confidence level, on the contribution to the signal 
from DM annihilation, XX v v, or decay, Y — > vv. We 
shall employ the DM decay limits of Refs. 2, 3]. A simple 
rescaling suffices to convert these limits into constraints 
on the decay mode of interest here, x ~^ x' v ^-i under the 
assumption that vv are emitted monoenergetically. The 
resulting constraints are shown in Fig. |3] 



III. LATE-DECAYING DARK MATTER AND 
THE SMALL-SCALE GALACTIC STRUCTURE 

The standard ACDM cosmological model successfully 
reproduces the observed structure of the Universe at large 
scales. However, at small scales, there appear to be 
several inconsistencies between observations and CDM 
simulations of galaxy formation. The negligible primor- 
dial velocity dispersion of CDM particles allows gravi- 
tational clustering to occur down to very small scales. 
Structure forms hierarchically, with small scales collaps- 
ing first and forming dense clumps. The resulting phase 
packing in the inner re gion s of galaxies leads to cuspy 
central density profiles 24], which are currently disfa- 
vored b y th e rotational curves of dwarf spheroidal galax- 
ies [iH I231 125rl28j . Dense lumps of matter also survive 
the mergers and become satellite galaxies. The num- 
ber of satellite galaxies predicted by CDM simulations 
greatly exceeds the number observed around the Milky 
Way [ii], H3, HH- A number of other disparities be- 
tween observations and CDM predictions have been re- 
ported. These include the overestimation of the number 
of halos in low-density voids [H, [35| , the non- prediction 
of pure-disk or disk-dominated galaxies |29| . the large 
angular- momentum loss by condensing gas [34j, and the 
"bottom- up" hierarchical formation |3l| . 

It is not yet known whether the solution to these in- 
consistencies lies outside the dark-matter sector. Indeed, 
individual astrophysical solutions to some of these prob- 
lems have been suggested. However, it is possible that 
the apparent discrepancies between CDM and observed 
galactic structure point toward a modification of the 
standard CDM scenario. The significance of these hints 
for deciphering the nature of DM is underscored by the 
fact that gravitational effects remain the only confirmed 
evidence for the existence of DM. 

A possible alternative to the CDM model is warm dark 
matter (WDM), whose thermal free-streaming proper- 



ties result in a suppression of structure on small scales. 
Although well-motivated particle-physics candidates for 
WDM exist, WDM lacks the "naturalness" of the stan- 
dard CDM scenario in which the observed relic DM abun- 
dance arises from thermal freeze-out (CDM candidates 
produced by other mechanisms also exist). It is possi- 
ble, though, to incorporate some WDM-like structure- 
formation features within the CDM paradigm, if we dis- 
card the notion that CDM particles are completely sta- 
ble. 

It has been suggested that if the heavy relic particles 
decay according to Eq. ([T]) , the energy released in the de- 
cay will change the standard picture of structure forma- 
tion. If the decays occur at early times, ~ 10 5 s — 10 8 s, 
before gravitational collapse, the massive decay prod- 
ucts acquire non-negligible velocities, which results in a 
suppression of the amplitude of density perturbations at 
small scales (still in the linear regime). The daughter 
DM particles behave effectively as WDM. Existing lim- 
its on the free-streaming and phase-packing properties 
of WDM can be directly translated to determine the in- 
teresting and the excluded regions of the CDM decay 
parameter space [HI, EH • 

Another interesting possibility, investigated in 
Refs. [IB-Gil, arises if the decays take place during 
the nonlinear stages of gravitational collapse. In this 
scenario the energy acquired by the decay products 
causes the halo to expand and, as a result, softens the 
steep central cusps predicted by CDM [16]. It also 
results in a decrease in the halo circular velocity, which 
can account for the deficit of satellite galaxies with 
velocity 10-20 km/s, with respect to standard CDM 
predictions [l7j . 

The efficacy and viability of the late DM decay mecha- 
nism to alleviate the CDM small-scale structure problems 
depend on two parameters: the parent particle lifetime, 
r, and the ratio e = Am/m x which gives the recoil ve- 
locity of the massive daughter particle, v x > ~ e. 6 Qual- 
itatively, one expects r to be comparable to the age of 
the Universe, so that the decays occur during galaxy for- 
mation. In the limit of very large r, the standard CDM 
scenario is recovered, while small r corresponds to the 
early decay scenario, examined in [Til Il5j. The velocity 
imparted to the heavy daughter particles should be suf- 
ficiently large - of the order of the virial velocity of the 
galaxy - in order to have an effect on the formation of 
the halo. It should not exceed, though, the escape veloc- 
ity from the halo, which would be catastrophic for the 
galaxy, unless, of course, r is very large. 

Abdelqader and Melia [lj| used a semi-analytical ap- 
proach that allowed them to incorporate the effect of 
DM decay in halo evolution. They showed that DM 



6 This becomes a strict equality in the case of a 2-body decay. If 
more than one relativistic particle is produced, as in the cases 
examined in this paper, the correction in v x i is insignificant. 
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FIG. 4: Constraints on DM lifetime r vs mass splitting Am 
and mass m x , for e = 10 -5 (recoil velocity v x / = 5km/s), for 
the decay channels x x' + 7 (g reen )> X X + e ~ + e+ 
(blue) and \ — > x' + v + v (purple). Color-shaded re- 
gions are excluded by the observed radiation backgrounds 
(Ref. Qj] and this work). The decays may have an observ- 
able effect on galactic halo structure if they occur at times 
0.1 Gyr < t < 100 Gyr (interval between solid lines). Longer 
lifetimes r > 100 Gyr correspond to the standard CDM sce- 
nario, while decays occurring at r < 0.1 Gyr (but after the 
matter-radiation equality) disrupt the galaxy formation sig- 
nificantly and are excluded [13, • 



decay preferentially heats smaller haloes, causing them 
to expand, and reduces their present-day circular ve- 
locity. They argued that if e ~ (5 — 7) x 10 -5 and 
r ~ (few — 30) Gyr, dark-matter decay can well account 
for the deficit in the observed number of galaxies with 
circular velocities in the range 10-20 km/s. 

Peter et al. [TH, [l9| performed a more extensive in- 
vestigation of the t — 6 parameter space, by means of 
both semi-analytical calculations and simulations. They 
demonstrated the way in which DM decay affects the halo 
mass-concentration relation and mass function, and used 
measurements of these quantities to constrain the DM 
decay parameter space. They located the allowed and 
potentially interesting parameter space roughly within 
the range 



e = 3.4 x 10~ 6 - 3.4 x 10~ 4 , 
r = (0.1 - 100) Gyr. 



(26) 
(27) 



The corresponding recoil velocity is v x > = (1 — 100) km/s. 
The actual limits on e and r are, of course, correlated. 

In Sec. |lll we derived constraints on the lifetime for 
late decay of DM into a massive daughter and either 
or a vv pair. For a given value of e, the bounds of Sec. [IT] 
determine whether these decays can take place in time to 
affect halo evolution. 

The flux of DM decay products in the Galaxy today is 
proportional to l/(mre T °/ T ), and the limits in Fig[T]arc 



FIG. 5: Same as in Fig. Q] but for e = 5 x 10~ 5 , or re- 
coil velocity v x i ~ 15 km/s, suggested in [17(] as suitable 
for resolving the missing-satellite problem. This would re- 
quire a dark-matter lifetime within the narrower interval 
1 Gyr < t < 30 Gyr shown (dashed lines). 
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FIG. 6: Same as in Fig. [4] but for s = 3 x 10 4 , or recoil 
velocity v x > ~ 90 km/s. 



expressed in terms of this quantity. Translating these 
limits to bounds on r alone, we see that the constraints 
may be satisfied for either (i) very long lifetimes, r > To, 
for which the current DM decay rate is very small; or 
(ii) very short lifetimes, t < rj, for which the abun- 
dance of the unstable parent in the Universe today is 
very small (as most decays have already taken place). 
Between these extremes limits lies a band of excluded 
lifetimes. These bands are shown in Figs. 0] - E] for 
three values of e spanning the range given in Eq. ([26]), 
e = 10~ 5 , 5 x 10~ 5 and 3 x 10~ 4 , corresponding to re- 
coil velocities v x i ~ 5km/s, 15km/s, and 90km/s respec- 
tively. 
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For the values of e considered, the excluded band spans 
several orders of magnitude in lifetime. These excluded 
bands become smaller as e is decreased, such that for suf- 
ficiently small epsilon there is no constraint. Note that 
while the upper limit of the excluded band may be im- 
proved with more sensitive flux observations, the same is 
not true for the lower limit, at t ~ 1 Gyr, since for these 
lifetimes the factor e~ T °/ T dominates. For this reason, 
the lower limit of the excluded regions for the photon 
and decay modes are very similar. 

We now compare the bounds on the DM lifetime r 
with the interesting range of values specified in Eq. (l27l) . 
Figures H] - [6] show that these bounds significantly con- 
strain, but do not rule out, unstable dark-matter models 
whose dominant decay channel produces the SM parti- 
cles considered. For the 7 or e ± decay modes, lifetimes 
greater than ~ 1 Gyr are eliminated, and we thus iden- 
tify t ~ (0.1 — 1) Gyr as the range of allowed lifetimes for 
which decays may affect structure formation. For decay 
modes to neutrinos, however, the allowed lifetimes span 
the entire interval of Eq. ([27]) (at least for some masses). 
Of particular interest, perhaps, is the case of weakly in- 
teracting massive particle (WIMP) dark matter of mass 
m x <~ 10 2 GeV, decaying via the vv channel, with r and 
e in accordance with the values suggested in Ref. [I?} 
for resolving the missing satellite problem. These decay 
parameters are indicated by the dashed lines in Fig. [5] 



IV. LIMITS ON MODELS WITH LATE 
DARK-MATTER DECAY 

A number of particle-physics models have been pro- 
posed which predict or invoke a DM decay mode de- 
scribed by Eq. (fTJ) [63 - l67j . The motivation for these 
models varies widely, thus the parameter space of inter- 
est does not always overlap with that relevant for struc- 
ture formation. We now describe, in general terms, when 
and how the bounds derived in Sec. [II] can constrain de- 
cay modes of the type in Eq. (JT]). We then apply these 
considerations to particular models found in recent liter- 
ature. 

In Secs.HTland lllll we assumed that all of the DM was in 
the form of the heavy, unstable, \ particles prior to decay. 
However, this need not be the case. In fact, many models 
predict or require that the relic \ abundance be only a 
fraction of the total CDM density. This may be true even 
if the decay time is larger than the age of the Universe, 
provided that either (i) other processes (e.g. inelastic 
scattering) efficiently eliminate the heavy particles (x) in 
favor of the light ones (%') [65[, or (ii) the x~ x' sector 
makes only a subdominant contribution to the total DM 
density. The relic x abundance is, of course, critical in 
deriving limits on specific models. 

Allowing the x abundance to be a free parameter, 
the four variables which determine the relevance of con- 
straints inferred from the present radiation backgrounds 
are (i) the mass splitting, Am, (ii) the mass of the DM 
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FIG. 7: The observed radiation backgrounds can constrain 
the lifetime of the heavy DM state only in models which lie 
in the regions above the lines. Models which lie below the 
lines cannot be constrained by the bounds presented in Sec. 
DTI As in Fig. [T] the different curves correspond to the decay 
channels x ~~ *• X + 7 (solid lines), x. ~~ * X + e ~ + e+ (dashed 
line), and x ~* X + v + v (dotted line). 



particle, m x , or equivalently e — Am/m x , (iii) the de- 
cay lifetime of the heavy DM state, t, and (iv) the relic 
fraction of DM in the form of x particles before decay 
occurs, f x = Y x /Yqbm- In terms of these parameters, 
the constraints of Sec. [II] can be expressed as 



Te To/T ^ {e- f x )T{Am)/Ar 



(28) 



where J(Am) can be read off Fig. [T] for the decay chan- 
nels considered. 



A. Fixed e ■ f x 

Consider constraints on the lifetime for models in 
which the parameters e and f x are specified. Similar 
to the analysis of Sec. IIIIl Eq. ([2"5]l implies an excluded 
band for r, provided that the fraction of energy released 
per decay (here quantified by s) and/or the fraction of 
decaying dark matter, f x , are sufficiently large. More 
specifically, in order for Eq. (|28p to translate into con- 
straints on t, the parameters e, f x must satisfy 



r e 



J(Am)/Am' 



(29) 



(since re T °/ T roe, for any r). 

In Fig. [7]we sketch the regions defined by Eq. (|2"9"]l . for 
the three decay channels under consideration. Models 
which lie above (below) the curves are constrained (not 
constrained) by the bounds derived in Sec. [II] using the 
observed radiation backgrounds. For models that lie on 
the curves, the constraint (l28l) produces a zero- width ex- 
clusion band for the lifetime at r = tq. Larger values of 
£■ f x yield exclusion bands which span a range of lifetimes 
around Tq. 
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The above discussion is relevant to a class of models 
proposed in Refs. [64], HH], in which DM possesses an ex- 
cited state. As a result of an approximate symmetry, 
the excited state is separated by a small mass splitting 
from the ground state, to which it decays with emission 
of pairs. For a particular implementation of these 
models, with Am « 2 MeV and m x ss 500 GeV, the ki- 
netic excitations and the subsequent decays of DM par- 
ticles inside the Galactic halo can produce the positrons 
which give rise to the observed 511 keV annihilation line 
in the center of our Galaxy 64]. For the range of pa- 
rameters favored by the 511 keV signal, f x < 10~ 3 and 
e ~ 10~ 7 — 10~ 6 , the scenario lies comfortably below the 
curve in Fig. and is thus not constrained by our 
observational bounds. However, significantly larger val- 
ues of e • f x are possible within the same generic class of 
models, for which our bounds would impose a non-trivial 
constraint. 



B. Fixed r 

Returning to the constraint of Eq. (|28p. we now con- 
sider models in which r (rather than e ■ f x ) is narrowly 
specified. Equation (|28fl then limits the energy release 
via decay, and relic abundance of the unstable state, as 
per 



J(Am)/Am' 
This constraint is substantial only if 

re ro/r < J 7 (Am) /Am 



(30) 



(31) 



(since by definition s,f x $J 1). We obtain a nontrivial 
constraint provided that 1 Gyr < r < 10 6 — 10 12 Gyr, 
where the exact range depends on the mass splitting and 
the decay channel considered. Models with r < 1 Gyr 
trivially satisfy the bounds placed by radiation back- 
grounds, as do models with very large r, in accordance 
with the discussion in Sec. IIIII 

The "degenerate gravitino" scenario of Ref. [67[ falls in 
the category of models whose energy release may be con- 
strained by observations. In this scenario, the lightest 
supersymmetric particle (LSP) and the next-to-lightest 
supersymmetric particle (NLSP) have a small mass split- 
ting, Am ~ 10~ 2 GeV, such that the slightly heavier 
NLSP can decay into the LSP with emission of a photon, 
with a lifetime comparable to the age of the Universe. 
These parameters satisfy the inequality (1311) and thus 
the observational bounds represent a non-trivial con- 
straint on the model. Specifically, the photon bounds [l| 
constrain the relic NLSP abundance before decay in a 
mixed (NLSP + LSP) DM scenario, as explicitly shown 
in Ref. 

As an example of models which trivially evade the 
constraint (|3T)]) . we mention here the inelastic DM sce- 
narios 66]. These models aim to reconcile the results 



of DAMA and CDMS direct detection experiments, by 
evoking small mass splittings Am ~ 100 keV. However, 
since the decays occur very early (before recombination) 
this scenario cannot be constrained by the present radia- 
tion backgrounds. Substituting the lifetime of the heavier 
state r ~ (10 2 — 10 3 ) yr into Eq. (|3~Tj) of course trivially 
verifies this fact. 



V. CONCLUSIONS 

The decay lifetime is a fundamental property of dark 
matter. Knowledge of this parameter may provide an es- 
sential clue for unravelling the particle nature of dark 
matter. In the standard CDM scenario, dark matter 
particle interactions play no role in the evolution of the 
Universe subsequent to DM freeze-out, which limits our 
ability to probe DM particle interactions via cosmological 
observations. However, there is a large class of DM mod- 
els in which DM particle interactions are in fact invoked 
to explain cosmological observations. One such hypoth- 
esis is that DM decay may provide a possible solution to 
the small scale structure problems of standard CDM. 

We have considered a class of DM models in which 
DM decays at late times to a nearly degenerate daughter 
plus relativistic particles, \ x' + I n this scenario, 
the relativistic daughters carry only a small fraction, e = 
Am/m x , of the parent DM mass and thus decays do 
not alter the DM energy density. However, the energy 
acquired by the daughter \' effectively heats the halo, 
and has been proposed as a means of modifying halo 
structure in such a way as to bring CDM predictions 
into better agreement with observation. 

We focused on the scenario in which the relativistic 
daughters consist of Standard Model particles, and de- 
rived constraints on the cases in which those daughter 
particles are either or neutrinos. In Fig. [T] we summa- 
rize these constraints and compare them with the anal- 
ogous photon bounds obtained in Ref. pj. Since pho- 
tons provide the most stringent bounds and neutrinos 
the least stringent, these results, taken together, delin- 
eate the range of lifetime bounds for these models. 

Our constraints on the decay mode lie between 
those for photons and neutrinos, though we note that 
there are some masses for which the constraints be- 
come comparable to those for photons. Note that many 
models that allow the process \ — > x'vv will also allow 
X ~^ x' e+e ~, f° r which stronger bounds apply. 

For the (and 7) decay modes, the lifetime must ei- 
ther be much greater than the age of the Universe (for 
which decay cannot affect structure formation) or smaller 
than ~ lGyr. Though this eliminates some of the param- 
eter space of interest in structure evolution, a significant 
portion remains open. For the vv decay modes, a much 
larger range of lifetimes is permitted. 

Models which are not necessarily motivated by struc- 
ture formation considerations, but nonetheless entail de- 
cay between particle states of similar mass, may also be 
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constrained by the bounds presented here. This is pos- 
sible, provided that the relic density of the heavier state 
from the early Universe is not too small, and/or the de- 
cay lifetime is larger than 1 Gyr. The constraints and the 
analysis presented here can provide generic useful guide- 
lines for constructing models which feature decay modes 
between closely degenerate states. 

We note that our limits in Figs. [T]- [3] are strictly appli- 
cable only for e <C 1 (so that the approximation Ei ~ Am 
is valid). However, for e = 1, our monoenergetic con- 
straints correspond to those for the decays x —> e + e~ and 
X — > vv (i.e., where there is no x')- For the vv case this 
simply reproduces the results of Ref. @j. For the e^- case, 
our bounds inferred using inverse Compton scattering 
and bremsstrahlung emission, synchrotron radiation, and 



the positron flux, reproduce the ones derived in Ref. [lOj , 
using the same response functions employed here. How- 
ever, the monoenergetic limits of Fig. [2] inferred from the 
in-flight and at-rest annihilations represent strong new 
constraints of general applicability for the decay mode 
X -> e+e~, in the interval 2 MeV < m DM < 1 GeV. 
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